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ABSTRACT 


A  radially  directed  load  is  suddenly  applied  to  a  portion  of  the  outer  surface 
of  a  circular  cylindrical  shell  which  responds  in  a  state  of  plane  strain.  An 
analytical  solution  for  the  resulting  dynamic  response  is  obtained  within  the  con¬ 
text  of  linear  elasticity  theory,  Flugge  shell  theox-y,  and  an  "improved"  shell 
theory.  A  comparison  of  the  analytical  solutions  ,.i.d  numerical  results  for  a 
specific  loading  indicate  that  the  improved  theory  is  far  superior  to  the  Fiiigge 
theory  in  terms  of  predicting'  both  the  magnitude  and  characteristics  of  the 
response.  However,  as  expected,  neither  shell  thooiy  can  satisfactorily 
predict  the  wave  character  of  the  initial  response-. 
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I.  INTRODUCTION 


The  small  motions  of  an  isotropic,  elastic  medium,  produced  by  a  dis¬ 
turbance  of  its  bounding  surfaces,  may  be  described  mathematically  by  the 
equations  of  the  thr<  e  dimensional  theory  of  elasticity  together  with  appropriate 
boundary  and  initial  conditions  [l].^  The  first  investigations  of  these  equations 
of  motion  in  cylindrical  coordinates  were  conducted  by  Pochhammer  [2]  in  1876 
and  Chree  [3]  in  1889.  Their  studies  dealt  with  the  propagation  of  free  harmonic 
waves  in  a  solid  cylinder  which  was  infinite  in  extent  in  the  direction  of  its 
generators.  Si  oe  then  several  extensions  and  refinements  of  these  initial  studies 
have  been  made,  most  notably  the  addition  of  numerical  data  for  the  frequency 
equations.  For  a  thorough  discussion  of  Pochhammer's  work  and  some  of  the 
subsequent  investigations,  the  books  by  Love  [  l]  and  Kolsky  [ 5]  should  be 
consulted . 

The  study  of  the  motions  of  cylindrical  shells  using  these  equations  is 
considerably  more  r<.^cnt.  It  was  only  in  the  past  two  decades  that  an  extensive 
effort  was  made  to  study  the  free  harmonic  vibrations  of  cylindrical  shells  as 
characterized  by  the  three  dimensional  theory  of  elasticity.  For  a  sample  of 
the  literature  on  this  subject  references  [6]  through  [19_]  should  be  consulted. 
Also,  recent  studies  have  been  made  of  the  forced  motion  and  transient  re¬ 
sponse  of  cylindrical  shells  using  this  theory.  vor  example,  in  1964  Liu  and 
Chang  [20l  investigated  the  transient  radial  displacement  of  an  infinitely 
long  cylindrical  shell  subjected  to  an  internal  axisymmetric  blast  load  and 

^ ^Numbers  in  brackets  designate  references  at  th<  end  of  the  paper. 


sudden  temperature  change.  By  applying  the  method  of  Mindlin  and  Goodman  [2l] 
to  the  problem,  they  "'ere  able  to  construct  a  solution  in  terms  of  the  noi.nal 
modes  of  vibration  of  the  cylinder. 

Subsequently  Suzuki  [22]  considered  the  problem  of  a  circular  ring  subjected 
to  a  transient  pressure  loading  of  both  the  inner  and  outer  surfaces.  He  attacked 
the  problem  with  a  combination  of  Laplace  Transforms  for  the  time  variable  to¬ 
gether  with  a  harmonic  analysis  with  respect  to  the  angular  coordinate  in  the 
plane  of  the  ring.  He  thus  formulated  the  general  problem  of  a  suddenly  applied, 
exponentially  decaying  load,  arbitrarily  distributed  over  the  lateral  surfaces  of 
the  ring.  However,  he  onlv  rented  solutions  for  the  axialiy  symmetric  case. 

In  1967  Garnet  and  Crouzet-Pascal  [23]  investigated  the  response  of  an 
infinite  cylindrical  shell  imbedded  in  an  infinite  elastic  medium  produced  by  a 
plane  dilatational  wave  traveling  through  the  medium  in  a  direction  normal  to  the 
cylinder's  axis.  Their  approach  was  to  construct  a  train  of  incident  pulses 
from  steady-state  components  such  that  each  pulse  contained  the  time  history 
of  the  transient  stress  in  the  incident  wave.  By  making  the  time  interval  between 
successive  pulses  sufficiently  large,  the  cylinder  would  return  to  its  original, 
unstrained  state  before  the  arrival  of  the  next  pulse  in  the  train.  This  occured 
because  of  the  radiation  of  energy’  from  the  cylinder  through  the  surrounding 
medium  to  infinity.  This  approach  proved  to  be  very  successful  and  results 
were  obtained  to  illustrate  the  time  history  of  the  stresses  and  displacements 
in  the  cylinder. 

The  examples  cited  above  illustrate  that  the  forced  motion  of  a  cylindrical 
shell  as  characterized  by  the  three  dimensional  theory  of  elasticity  is  mathe¬ 
matically  very  complex.  A  quantitative  description  of  the  response  is  extremely 
difficult  without  the  use  of  high  speed  computers.  This  is  the  principal  reason 
for  the  long  delay  between  the  initial  investigations  of  Pochhammer  and  Chree 
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and  those  just  mentioned.  The  cylindrical  shell,  however,  is  *  very  common 
element  with  many  and  varied  applications,  therefore  the  need  to  analyze  its 
dynamic  response  arose  long  before  the  means  for  carrying  out  such  an  analysis 
within  the  framework  of  the  three  dimensional  theory  of  elasticity  were  avail¬ 
able.  This  lead  to  the  development  of  several,  mathematically  simpler, 
theories  to  describe  the  motion  of  cylindrical  shells.  These  so  called  shell 
theories  were  based  on  the  assumption  that  the  radial  thickness  of  the  shell 
was  much  smaller  than  the  radius  of  the  nedfan  surface  of  the  shell.  With 
this  assumption  the  dependent  variables  could  be  3xpanded  into  convergent 
power  series  in  the  thickness  coordinate  and  the  first  one  or  two  terms  in  these 
expansions  would  suffice  to  describe  the  response  of  the  shell.  A  theory  of  this 
type  was  developed  by  Love  [4]  at  about  the  same  time  as  the  Pochhammer  and 
Chree  investigations.  Since  then  numerous  other  shell  theories  have  been 
proposed.  Most  of  these  may  be  placed  into  one  of  the  following  three  categories. 
The  first  type  of  theory  is  called  a  membrane  theory.  Here,  n  >  variation  of 
the  dependent  variables  through  the  thickness  of  the  shell  is  permitted.  Re¬ 
ference  to  this  type  of  theory  is  made  by  Rayleigh  [24]  and  the  equations  of 
motion  for  a  cylindrical  shell  may  be  found  in  the  books  by  Flugge  [25]  or 
Vlasov  [20] .  The  second  category  contains  the  classical  shell  theories.  These 
allow  the  dependent  variables  to  vary  linearly  through  the  thickness  of  the 
shell  but  in  such  a  manner  that  straight  line  elements  normal  to  the  median 
surface  of  the  shell  in  the  unstrained  state  remain  normal  during  the  motion 
of  the  shell.  Furthermore,  these  elements  retain  their  original  length  and 
contribute  no  rotatory  inertia  to  the  motion.  Love's  equations  are  contained 
in  this  category  along  with  those  of  Flifgge,  Donnell,  Vlasov  and  Sanders. 

The  third  category  contains  the  improved  theories.  As  in  the  classical  theory 
the  dependent  variables  are  allowed  to  vary  linearly  through  the  thickness  of 
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the  shell.  The  improvement  is  obtained  by  allowing  the  previously  mentioned 
line  elements  to  rotate  relative  to  the  median  surface  and  by  including  their 
rotatory  inertia  in  the  motion  of  the  shell.  Although  these  line  elements  are 
still  required  to  remain  straight  and  retain  their  original  length,  a  further 
improvement  is  obtained  b\  introducing  a  correction  factor  into  the  trans¬ 
verse  shear  force  to  compensate  for  this.  The  magnitude  of  the  correction 
factor  is  obtained  by  matching  the  phase  velocity  of  the  lowest  mode  of 
propagation  of  frci  harmonic  waves  in  the  axial  direction  with  that  obtained 
from  the  three dimensional  theory  of  elasticity  11 143  .  The  equations  of  motion 
of  a  cylindrical  shell  characterized  by  the  improved  theory  may  be  found  in 
the  papers  by  Herrmann  and  Mirsky  [ 27  j  and  also  Reismann  and  Medige  [28]. 

As  a  result  of  the  diversity  of  the  proposed  shell  theories,  the  following 
question  arises.  For  a  cylindrical  shell  subjected  to  a  specific  disturbance, 
which  of  these  theories  predicts  the  response  to  within  a  given  accuracy  with 
the  least  effort0  This  question  may  be  answered  by  comparing  the  response 
predicted  bv  each  of  the  shell  theories  to  the  response  predicted  by  the  three 
dimensional  theory  of  elasticity  for  each  specific  disturbance.  This  approach, 

however,  would  negate  the  only  advantage  of  the  shell  theories,  which  is  their 
relative  simplicity  compared  to  the  three  dimensional  theory  of  elasticity.  An 
alternative  approach  is  to  carry  out  a  comparison  with  the  elasticity  theory  in 
only  a  few  specific  eases  which  represent  the  limits  of  the  range  of  possible 
disturbances  and  shell  geometries.  From  these  few  limiting  cases  rational 
estimates  oi  the  accuracy  of  the  shell  theories  for  various  other  disturbances 
and  shell  geometries  could  then  be  made.  Some  of  these  limiting  cases  have 
already  been  investigated  and  comparisons  of  the  shell  theories  with  elasticity 
theory  have  been  made.  For  example  Klosner  in  references  f  20  j  through  [32] 
and  Iwngor  and  Yogananda  .33’  have  examined  various  problems  involving 
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cylindrical  shells  which  are  statically  loaded  using  elasticity  theory  and  shell 
theories.  However,  no  comparisons  have  as  yet  appeared  for  the  other  extreme 
ease,  that  is  where  the  load  is  suddenly  applied.  Since  most  real  loading  situations 
fall  somewhere  between  these  two  extremes  a  comparison  of  the  shell  theories 
to  elasticity  theory  for  the  latter  case  would  be  of  great  interest.  Therefore,  the 
primary  purpose  of  this  investigation  will  be  to  present  such  a  comparison. 

To  accomplish  this  goal  a  cylindrical  shell  of  circular  cross  section  which 
is  infinite  in  the  direction  of  its  generators  is  subjected  to  a  suddenly  applied 
force  on  its  outer  surface.  This  force  is  chosen  to  act  only  in  the  radial 
direction  and  also  to  be  invariant  in  the  direction  of  the  generators  of  the 
cylinder.  With  these  restrictions  on  the  force  the  cylinder  may  be  assumed  to 
be  in  a  state  of  plane  strain  and  therefore  only  a  plane  section  normal  to  its 
generators  will  be  considered.  The  response  of  the  shell  predicted  by  the  three 
dimensional  theorv  of  elasticity  will  be  found  using  the  method  discussed  by 
Suzuki  i.22l.  Next,  the  response  of  the  shell  predicted  by  the  Fliigge  Theorv 
f  25 a  classical  theon  ,  and  also  the  improved  theoi’y  due  to  Herrmann  and 
Mirskv  [27]  will  be  found  from  their  corresponding  Green's  functions  given 
by  Pawlik  and  Reismann  in  [  ,34 1. 

The  time  history  of  the  displacements  and  stresses  at  various  points  in 
the  shell  as  characterized  by  the  throe  theories  will  then  be  compared  for  a 
specific  shell  geometry  and  load  distribution. 
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An  infinitely  long,  circular  cylindrical  shell,  composed  of  an  isotropic, 
linearly  elastic  material,  is  at  rest  and  in  an  unstrained  state.  A  load,  directed 
radially  inward,  is  suddenly  applied  to  the  outer  cylindrical  surface  of  the  shell. 
The  distribution  of  the  load  on  this  surface  is  constant  along  the  generators  of  the 
cylinder  but  otherwise  arbitrary. 

To  simplify  the  analysis  of  the  resulting  motion  of  the  cylinder,  its  weight 
will  be  neglected.  In  this  case  it  may  be  assumed  that  the  motion  of  any  plane 
section  normal  to  the  generators  of  the  cylinder  takes  place  entirely  within  thai 
plane  and  is  identical  for  ill  such  sections.  In  other  words  the  shell  responds  in 
a  state  of  plane  strain. 

The  initial  geometry  of  a  plane  section  of  the  shell  together  with  a  particular 
load  is  shown  in  Figure  1 .  The  applied  load  may  be  represented  symbolically  as 
follows . 

P(0,t)  -  PQ  g(@)  H(t)  (1) 

In  (1),  PQ  is  a  reference  pressure,  g(0)  describes  the  distribution  of  the  load  over 
the  outer  boundary  of  the  shell  and  H(t)  is  the  Heaviside  step  function  with  respect 
to  time: 


H(t) 


0  ,t'0 
1  .  t>  0 


(2) 


The  analysis  of  the  motion  for  an  arbitrary  load  distribution,  g(0),  is 


facilitated  by  expressing  g(9 )  in  terms  of  a  Fourier  series. 


1  *() 

g(6)  ~  j_  ~  '  (a^cos  n0  -  b^sinnS) 

n  1 


rn 


t 


However,  since  ultimately  a  quantitative  measure  of  the  response  is  desired,  the 
above  series  will  necessarily  be  truncated  after  summing  some  finite  number  of 


terms  (N).  This  results  in  an  approximation  of  the  actual  distribution  g(9)  by  the 


finite  sum  gN(0).  The  accuracy  of  this  approximation  for  any  given  value  of  N 

may  be  increased  by  applying  the  process  of  "smoothing, "  explained  in  Appendix 

to  fN(P)*  In  the  present  case  this  amounts  to  multiplying  each  term  of  the 

truncated  series  by  the  factor  ( —  sin  ~  ). 

vnrr  N  / 


(ancos  nS 


n  =  1 


(4) 


where  a  =  a  sin  ~ ) 
n  n  xmT  N  / 

a  =  b  (X  8ln  m) 

n  n  snn  N  / 

The  analysis  will  therefore  be  carried  out  for  the  load  distribution  gN(0) 
given  by  Equation  (4), 

All  the  variables  used  in  the  following  analysis  are  in  dimensionless  form 
both  for  convenience  and  generality.  The  conversion  to  dimensional  form  is 
given  in  the  Nomenclature. 
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IH.  THEORY  OF  ELASTICITY 


A.  Basic  Equations 

Navier's  equations  of  motion  for  an  isotropic,  linearly  elastic  material 
may  be  written  in  dimensionless  form  as  follows. 


y2v2 u  +  (l  -y2)$  <?-u)  = 

at2 


(5) 


Here  u  is  the  displacement  vector,  V  is  the  gradient  operator,  v2  =  V  •  V  is 
the  Laplician  operator  andy  is  the  wave  speed  ratio.  Using  Helmholtz's 
Theorem  [36]  the  displacement  vector  field  may  be  expressed  as  follows. 

U=^0+Vxt,V»f*=O  (6) 

<0  is  the  scalar  potential  and  iji  the  vector  potential  of  the  vector  field  u .  Sub¬ 
stituting  (6)  into  (5)  results  in  the  following  equation. 


#  [A '  +  ? x  jy2?2 $  -  =  0 

at  at 

This  equation  is  satisfied  if 

ar 


and 


A2?.*- 

ar 


(7) 


For  the  case  of  plane  strain  in  cylindrical  coordinates  the  displacement  vector 
and  potentials  are 

u  w(r,8,t)  sr  ♦  v(r,  8,t)eQ 

o  -o(r,0,t)  (8) 

1  ;  (r,  0,t)  e 
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where  e  ,  $  are  the  unit  vectors  in  the  radial,  circumferential  and  axial 

v  o  z 

directions  respectively.  The  kinematic  and  constitutive  relations  for  the  plane 
strain  case  are 
dW 

er  "  >r 

1  3v  w  /n. 

e0  =  r  }  0  4  r  (9) 


1  J_W  t 

r  3  0  ‘  3r  r 


'r  "  er  ‘  (1  ~  2y>8 


0  "e0 


(1  -  2y2)( 


T  -  y  s 

where  e  .  ea,  -  ,  i.  are  the  normal  strains  and  corresponding  stresses  while 
r  a  r  y 

8  and  t  are  the  shear  an^le  and  shear  stress  respectively. 

Substitution  of  (8)  into  ((>)  and  (7)  yields  the  following  set  of  equations: 


r~ 

i 

♦  — 

i 

,2 

i 

0 

r- — i 

L  -  2 

r 

ir 

2 

r 

i02 

dt2 

2  r  *2 

.  i 

1 

2 

'1 

-  ]  ,  , 

.2 
o _ 

i  0 

r 

hr 

2 

r 

902 

J  * 

At2 

O  1  '' . 

W  1  -  —S 

t  r  r  >0 

(12) 

V  -  "  "ft  - 

r  i0  -ir 

Substituting  (12)  into  (9),  then  (S»>  into  (10)  results,  with  the  aid  of  (li),  in  the 
following  relations: 

„  220„2y2fl_^)  J_  A**  _ 

r  .  .2  “y  Lr  ir  2  n2  r-  Vr  tWJj 


„  JLO_2  2f  IJs2.t  -L  JL£  _ 
r  dt2  “y  Lr  ir  f2  nfl2  r- \v  W Jj 

,  _L2s^  .  9V2fiLP  ,  2_  f  1  JLL  jl 

8  Ot2  kr2 
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T 


(13) 


a‘.  *trA  a  /'I  w\1 

n2  ‘2y  LSr2  arlr  39  U 


The  initial  conditions  for  the  problem  under  consideration  are 
<0(r,9,O)  =  l(r,9,0)^  0 


d<3 


t  =  0 


=  0 


The  boundary  conditions  are 
<Tr(X,0,t)  -  -  gN(0)H(t) 
r(X,9,t)  -  0 
rrr(Y,0,t)  =  O 
r(Y,0,t)  =  0 


(14) 


(15) 


where  X  =  x  =  1  +  —  and  Y  =  y  =  1  -  are  the  outer  and  inner  radii  of  the  shell, 
respectively. 
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B.  Solution  of  the  Equations  of  Elasticity 

The  solution  of  the  system  of  equations,  boundary  conditions  and  initial 
conditions  (11)  through  (15)  will  now  be  obtained  by  a  combination  of  harmonic 
analysis  and  Laplace  transformation. 

In  view  of  the  boundary  conditions  (15)  and  the  form  of  g^T(0)  given  by  (4), 
the  dependent  variables  in  (11)  through  (15)  may  be  assumed  to  be  in  the  following 


form. 


ra  1 

o(r,0,t)  -  --  Lt  '  L  $n(r»t)(ar"os  "9  “0nsin  n0)J  <16_a) 


:  (r,e,c) 


'J'n(r>t)(ansin  n0  *  0ncos  n0) 


(10-b) 


i  ra 

i  r  o 


w(r,0.t)  *  ~  L'2  Wn(r,t)  •  _  Wn(r, t)(ancos  n0  -  ^sin  n9)J  (1?_a) 

n  1 


v(r,0,t) 


V  (r , t)(.r  sin  n9  fi  cos  n9) 
_  n  '  n  n  ' 


(17-b) 


-r(r,e,t)  -4  [_-xso(r‘t) 


r  1 

S  (r,t)(a  cos  n 0-/9  sin  n9)  (ls-a) 
n  n  Kn  J 


— n  — 

)(r,9 ,  t)  -4[  ~T  So(r’t>  ‘  ^  S®  (r,t)(.incos  n0  -  sin  n9)  1  <l>-b) 


(r,0,t)  ^Tn(r,t)(insin  n9  •  0ncos  n9) 


(lH-c) 


The  modal  coefficients  3^,  etc.;  are  determined  by  substituting  the 
assumed  solutions  (Hi)  through  (Is)  into  the  relations  (11)  through  (l.">).;  This 
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results  in  the  following  set  of  equations. 

.  11  n!l4  !\ 

r  Sr  r2_|V  a,2 

rsr  r2j\  y2  S(2 

%(r,t)  =s  0 
*n(r,0)  =  *n(r,0)  =0 

_ H _ n 


sJ(X,t)  *  H(t)  ,  sJ(Y,t)  =  0 


(19-a) 

(lJ-b) 

(19-c) 

(20-a) 

(2C-b, 

(2l-a) 


Tn(X,t)  =  0  ,  Tn(Y,t)  =  0 


wn<r-> 


n  n  . 

— —  t  --  * 

dr  r  n 


V  (r,t)  -  — -2-  t  ^<|> 
n'  ’  dr  r  n 

sr.r  ..  *\  2*f 


r  ( Wn  '  nV„] 


2 

Tn(r’t>‘^T--i?-Lv„-nWn] 

g  9  d  "<t»  0  2 

^r.t)  -  (1  -  2y“)  — ~  -  -sj-  [W  -  nV  ] 
11  dt  r  n  n 


where  n  =  [0,1,2,  ...  ,N}. 


<21— b) 
(22-a) 
(22— b) 
(22-c) 

(22— cl) 
(22-o) 


The  solution  of  the  differential  equations  (19)  may  be  obtained  by  application 
of  the  Laplace  transformation  defined  below. 

T(r,P) 

o 
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(23-b) 


f(r,t)  --£i  j  f(r,  P)e  ld. 

C-i‘° 

In  the  above  integral,  known  as  Bromwich's  integral  formula,  i  -  J-i  and  the  real 
number  C  is  chosen  so  that  P  -  C  lies  to  the  right  of  all  the  singularities  of  f(r,P) 
in  the  complex  P-plane. 


Applying  the  transformation  (23-a)  to  the  set  of  equations  (19)  through  (22) 


results  in  the  following  system  of  equations: 


=  0 


(24-a) 

<21-b) 


fo(r,P)^ 

0 

(24-c) 

SJ(X,P)  - 

1 

P  • 

^(Y,P)  =  0 

(25-a) 

TnvX,P)  = 

'J  , 

T  r  =  o 

ri 

(25-b) 

d$ 

n  •s’ 

1  r\ 

Wn(r,P)  ^ 

n 

dr 

(26-a) 

d* 

, 

r  n 

V*’P> 

_ n_ 

dr 

(2fi-b) 

^(r,P)  - 

2  T 

P  $  - 

2''2  -  , 

•  "J- -  Wn  '  "V 

(£li-C) 

Tn(r,P) 

p2$ 

n 

9V2 

-^Vn-nWn] 

(20-  ) 

3^r’p) 

(1  -  2y 

V2$„  -  -^-[W  -  -VI 

n  r  n  n 

(26-e) 

where  n  -  { 0, 1, 

2,  .  .  . 

.  n)  . 

The  solution  of  equations  (2-4)  is  (see  for  example  [35]) 

%(r,P) 

A  (P)I  (Pr)  *  B  (P)K  (Pr) 
xv  '  n  nv  '  nv 

(27-a) 

*n(r,P) 

C  (P)I 
n 

i  ~  j  -  D  (P)K  l— ) 
n  \  y  /  tv  ’  n\y  J 

(27-b) 
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(27-c) 


^0(r,P)s  o 

where  I  and  K  are  the  n  'dified  Bessel  functions  of  the  first  and  second  kind 
n  n 

respectively.  The  constants  A.  ,  B  ,Cn  and  D  are  determined  as  follows.  Sub- 

n  n  u  n 

stitute  (27)  into  the  relations  (26)  to  obtain  expressions  for  the  modal  stresses.  Nexi, 

substitute  these  expressions  into  the  boundary  conditions  (25).  This  results  in  a 

set  of  linear  algebraic  equations  from  which  the  constants  A  ,  B  ,  C  and  D  may 

n  n  n  n 

be  uniquely  determined.  After  substituting  the  constants  thus  obtained  into  (27)  it 
io  observed  that  all  the  modified  Bessel  functions  may  be  conveniently  grouped  into 
four  new  functions  called  the  cross,  products  o?  the  modified  Bessel  functions.  These 
are  defined  below: 


F^P.X.Y)  3  In(PX)K. (<PY)  -  In(PY)Kn(PX) 
F^P.X.Yt  ?  PYtln(PX)K^(PY)  -  I^(PY)Kn(PX)] 
F^P.X.Y)  =-  PX[l^(PX)Kn(PY)  -  yPY)K'(PX)] 
l'(PX)K'(PY)  -  I'(PY)K'(PX) 


fJ^P.X.Y) 


,  n  -  0 


P2XY[I^(PX)K^(PY)  -  l'(PY)K'(PX)]  ,  n^  0 


(28-a) 
(28— b) 
(28-c) 

(28-d) 


where  a  prime  denotes  differentiation  with  respect  to  the  argument  of  the  function. 
The  properties  of  these  functions  used  in  the  forthcoming  analysis  are  listed  in 
Appendix  II. 

In  terms  of  these  cross  products  the  solution  may  now  be  written  as 


Cftr.P) 

$(r,P)  7T~ - 

PJDn(P) 

Cn(r,P) 

Vr*P)  •  — 

p3n  (P) 


>17  0 

o 


(29-a) 


(29-b) 


where  C^anu  C‘  are  given  below, 
n  n 

For  n  -  0: 
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(30-a) 


2 

5j(r,P)  -  p^F^(P,r,  Y)  -  f£2)(P. r,Y) 

2 

Dq(P)  -  P2F^)(P,X,  Y)  -  F^P.X.Y) 


-  F^P.X.Y)  <  -^-F^P.X.Y)  (30-b) 

y 

For  n  =  1:. 

Cf(r,P)-  —  A^r.P)  (31-a) 

i  p~  i 

cJ'(r.P)  *  --A^r.P)  (31-b) 

DX(P)  -  (P)  (31  -c) 

For  n  -  {2,3,4 . N} : 

g£(r,P)  A(^(r,P)  Mn2- l)B^(r,P)  (32-a) 

£„(r,P)  -  Sn(r,P)  4  (n2-l)B';(r,P)  (32-b) 

Dn(P)  -  A*(P)  *  (n“-  )B*(P)  (32-c) 


—  o 

The  functions  A^  etc.  are  defined  below  in  terms  of  the  cross  products  for 
n  =  {1,2,3,  .  .  .  ,  Nj . 


a>) = *44-  ■  [P2  • 

n  x“V“  X  X^  ; 


Y  2  J  n 


-  ^l1’2  '  i^r_^n'2,(p)  -  1  J^F-Ai1'3><P» 

^{fn  ,4><P)  ‘  Fif'3)<P>] 


X  Y 


B^(P) 
n  ' 


1)  .  2  22  ,  1  1 
— 1  ■  4y  n  P  l -s-  * 

Y  Xw 


4y^  ( fry  *n2(n2- 

pW  \  x2y2 

■  n2^^2^  .  4V2„2P2  A  ■  V;  ■ '  pl  V4 ' 4  !><’■  "0-1 
L  x-y-  V  X"  V  X'J 


(continued  on  next  pane) 
l(i 


&  £3 


where  F 


J. 


X2r4v4n2ln2-  1)  .  4y2nV  +  p4lf(2,2)  p 

v4  Y2  J  " 


X 

r2  I-4-.4_2._2 


Y2  f 4vV(n  - 1)  4y  n  P  p4lp(3,3).p) 
‘  y4  Y2  J" 


+  4y4(n2-..4)..  f(4'4)(P)  -  ■4y--!^—  F(1,2)(P) 
.*2. .2  n  11 


x2y2 


-  4y2n2p2-Fj11'3>(P)  *  4y2P2[-ijF<n3'4)(p)  *  ^f[12'4)(P)J 


V-2 

Xs 

V-2 

y2 

y-2 

Xs 


Y2 


,*•%>  •  I[f<>.X.Y,F«>(2  X.y)+  fj>(£. X.Y^«>(P,X.Y)J 
[r,P) »  -  ^V^V.X.r)  -  F<n3)(P,X,r)] 

4  (p2  +  )  (p2  4  J*¥-)Kl><p-r- ’(? ,X-y) 


2*1 1 
Y2 


(p2  +  ^4^-)[F|i1|(P,r,Y)Fji2»tf,X,Yj 
X 


*  f<2,(P.F.V)F<1)(f,X.Y)].^-(p2+  ^)F'I»«P.r.Y,F<3>(i,X,Yj 
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*  ^K>(E  ^)nV.X,Y)  ♦  F<2>(f  .r,Y)p<>,X,Y,J 

5n(^^^fA2[r«3>(P,X,r,-F«)(P,X,)j 

*  n2(P2X2  ^  2y2n2)TW(P,r,Y)F^(|,X,Y) 

-  2y2n2F^(P,r,Y)F^(-E,X,Y)  ♦  2y2F®(P,r,  Y)fW(|,X,  y) 

-  (P2X2  4  2y2„V|i2>(P,r,Y)F<n2>(|,X,Y)| 

{2F2n3>(f  .M-  (A2  4  2y2„2)F<‘)(|,X,r) 

'  2F2”2ri‘\f.F,Y)FW(P,X,Y)  -  2y2n2F<1)(|,r,Y)F®(P,X,Y) 

-  2F2P'2\v.r.Y)Ff(P,X,Y)  *  2y2lf  \f  .r.Yj^V.X, Y)| 

The  above  relations  will  be  referred  to  as  equations  (33). 

The  solution  in  the  time  domain  will  now  be  obtained  by  substituting  (29) 
into  the  Bromwich  integral  formula  (23-b). 

....  ,  fCtl"  C?(r.P, 


V’**  YTi 


Vr-‘> 


0  Pt1T1 

— o - e  dP 

P  D  (P) 

>  n'  ' 

c  i  i  «  -  ;• 

C",r'P)  «HP 

~o - e  dP 

P  D  (P) 
n'  ' 


(34-a) 


(34-b) 


For  r  >  0,  each  oi  the  functions  C£(r,P),  ^(r.P)  and  Dn(P)  is  (a)  an  entire 

function  of  P,  (b)  symmetrical  in  P  and  (c)  nonzero  at  P  0.  D  (P)  has  a  do- 

n 

numerable  infinity  of  simple  zeroes  located  along  the  imaginary  axis  of  the  complex 

P -plane  for  each  value  of  n.  If  the  magnitude  of  these  zeroes  is  denoted  by  <*>  then 
_  -  nj 

D  ft  i  to  ,)  0  ;  i  0  12  /or  v 

n'  ,\y  ’  •'  (35-a) 


IK 


0 <  co  <  oo  ,  <{*)„<»•* 
no  nl  n2 


(35-b) 


The  integrals  (34-a,b)  will  now  be  evaluated  by  contour  integration.  From 
the  above  discussion  we  conclude  that  both  integrands  have  simple  poles  at 
P  =  ±  i  oo nj  and  a  pole  of  order  three  at  the  origin.  By  applying  the  residue  theorem 
to  the  integral  around  a  Bromwich  contour  with  C>  0  the  solutions  are  obtained 
in  the  following  form. 

For  t  >  0: 


*„<*■» -«£V.t>*  I 

j  -0 

33 

*n<r,t).*<f>(r,t)4  I  C*(r,Un))Qn)(t, 

j  =0 


where 


2  c® 


<’ . 

n'  ' 

(Si  if  a2  /Cn<r’P)ePtM 

n  ZLdP  '  D  (P)  /JP  =  ( 

n'  ' 


(36-a) 


(36-b) 


(37-a) 


(37— b) 


c£(r,oo)  -  C^(r,  ±  ico) 


(38-a) 


CWr.ttj)  -  C;(r,  i  iu>) 


(38-b) 


Dn(w)  "  Dn(±  iW) 


Q  (t) 

nj'  ' 


cos  oo  t 
nj 


oo  ,G  (oo  .) 
nj  n'  nj' 


(38— c) 


(39— a) 


Gn(W)  Gr(±  ico) 


5„(P) 


,  JD  (I>) 
1  nv  ’ 

2P  dP 


(39-b) 


(39-c) 


Bv  substituting  the  potentials  (3(>)  into  f"l)  we  obtain  the  desired  solutions  for  the 


mod';!  displacements  and  stresses. 
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W  (r,t)  =  W(s)(r.t)  +  V  W  .(r)Q  .(t) 
n'  '  n  '  ’  '  u  nj'  '  nj' 

j  =0 
00 

l  Vnj(r)Qnj(.) 

j  -0 


S>-‘>  =  Sn  9)<r’*>  *  2.  SnVr)Qn|<t> 

j  =0 


CO 


sJ(r,t).S®'s)(r,t)+  2  4,r)V) 

j-0 


■yr.y.T^W).  2,  TnJ(r)Qn.(t) 

j  =0 


CO 

r 


where 


W  Ar)  -  ~-C°(r,(ji  .)  +  -C*(r,  ta)  ,) 

nj'  '  3r  n'  nj7  r  n'  *  nj; 

V  ,(r)  =  .)  +  -C^(r,  u)  ,) 

nj'  '  dr  n'  nj'  r  n'  nj' 

2 

Sr,(r)  =  -  oo2.C0(r,oJ  .)  -  -^-[W  .  -  nV  .] 
nj'  1  nj  n'  *  ny  r  L  nj  nj 

2 

S®  (r)  -  u2.(l  -  2y2)C!V.u>  .)  +  -^-[W  .  -  nV  .] 

nj  '  ny  '  ’  n'  nj'  r  nj  nj 

2 

T  .(r)  -  -  uj2.C  (r,ou  )  -  ^ — [V  .  -  nW  .] 
nj'  ’  n.)  n'  ’  nj'  r  ^  nj  nj 


(40-a) 


(40-b) 


(40-c) 


(40-d) 


(40-s) 


(41-a) 
(41-b) 
(41 -c) 
(41-d) 
(41-e) 


O  lf 

The  functions  C^fr.ou),  (r ,  oo)  etc.  may  be  conveniently  expressed  in  terms 
of  the  cross  products  of  the  Bessel  functions  defined  below. 


Fjj  '(uj.x.y)  -  Jn(oux)Yn(4>y)  -  Jn(joy)Yn(u)x) 
F^Wx.y)  =  u!\’r Jn(olx)Y'(u)y)  -  4'(^v)Yn(jOX)] 
F^fcu.x.v)  -  jux'.l'(u>x)Yn(ujy)  -  -ywvjY'fujx)’ 


n 


VoU*)\'0U»  -  ^(wvlY^x)  n  0 

ufx\  r  j'(u;x)Y'(u;v)  -  .l'(a;\)Y,(jux)'  n  0 
n  n  n  n 


(42-a) 
(42-b) 
( 42-c> 

(42-d) 
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where  a  prime  denotes  differentiation  with  respect  to  the  argument  of  the  function. 
The  properties  of  these  functions  used  in  this  analysis  are  given  in  Appendix  II. 

The  functions  appearing  in  the  modal  solutions  (40)  and  (41)  are  given  below 
explicitly  in  terms  of  the  cross  products  of  the  Bessel  functions  for  the  three 
cases  n  =  0,  n  =  1  and  n  2  2. 

For  n  -  0  the  eigenvalues  (ooQ^;  j  -  0,1,2,  ... )  are  the  positive  real  roots 
of  the  equation 

D0M  -  0.  (43) 

The  functions  used  in  the  modal  solutions  are 

Do<«»  =y[w2FoW.y>  4  F^fco.x.y) 

y 

+  ^VFo3)(uJ*x,y)  4  (44-a) 

X 

2 

Cq(i\u))  -  ^  u)2F^(<d,r,y)  +  -^-F^fui.r.y)  I  (44-b) 

y  -J 

rC^fr,  UJ)  =  F^Wr.y)  •  -^FJ>4)(a),r,y)j  (44-c) 

G0(4  -i[(l-2y2)F^(cu,x,y)  *  i^l  -fl)F^(u;,x,y) 


*  K1 ‘4V)fo3)(^x>v) 

U!  X 

1  /.i  VV  •» 


(44-d) 


W^r.t) 


sfV.t) 


2>2(i-y2)d  -  V) 

w 

X 

2 

1  - 

2 

r 

2 

1  - 

1  o 


(15-a) 


(45-1)) 
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(45-C) 


2 

i*  h 

S®(s>(r,l)  — V 

*-  V 

x 

For  n  =  1  the  eigenvalues  (u)j ^ ;  j  =  0, 1 , 2,  .  .  . )  are  the  positive  real  roots 
of  the  equation 

D^(co)  =  0.  (46) 

The  corresponding  functions  used  in  the  modal  solutions  are 


D^co)  =  -  -VAf(W) 


-1_  A^/ 

2  Al' 

Hi 

J_A>| 
2  Ar 

Hi 


9  0® 

1 

dr 

2  1 

Hi 

9C* 

1 

A  t  / 

dr 

2  "1 

Hi 

Gi<“>  *  ^[Af  ,m)  -  iAN] 


(47-a) 

(47-b) 

(47-c) 

(47-d) 


where  a  prime  denotes  partial  differentiation  with  respect  to  r  and  a  dot  denotes 
differentiation  with  respect  to  hi-  A^5,  A’  ,  A^etc.  will  be  defined  shortly. 


vfV.t)  •  fj- 


2  2  , ,  „  2  <  2  N 

4  _ * _ i X-  .  ( \J—) 

2  2  \  2  2  /  2  / 
x  +  y  r  1  -  y  y 

ay  ) 

l  y  y 

2  2  „  2,2. 

_  _* _  (X-  ,  i  -L_  ) 

2  2  \  2  L  2  /  2  7 

x  •  y  r  1  -  y  y 


(48-a) 


In  -  -  2 


c-/-;-  3j 


(48-b) 
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x  +y  r  y /J 

(48-c) 

X  +  y  >  J  *  m* 

(48-d) 

T<s>(r,t)=4S^r(r--!!)+  x2  ixl.ill 

1  4a  L\y  r  J  x2+y2  V  r3  y  JJ 

(48-e) 

For  n  =  [2,3,4 . N]  the  eigenvalues  j  =  0,1,2,  .  . 

. )  are  the 

positive,  real  roots  oi  the  equation 

Dn(us)  -  0. 

(49) 

The  functions  needed  for  the  modal  solution  are 

D„<*>  A*(w)  +  (n2-  l)B^(uj) 

(50-a) 

C^(r,W)  -  A^(r,<4  -  (n2-l)B^(r,u)) 

(50-b) 

Cn(r,w)  "  An<r*W)  +  (n2-l)B*(r,w) 

(50-c) 

bC° 

(50-d) 

dCn 

dr  =  A,i#(r.w)  '  (n2-l)Bn'(r,<4 

(50-e) 

°n(w)  *  ‘  2^  ^An  (c0)  +  <n2-1)®^(^ 

(SO— f) 

where  a  prime  is  used  to  denote  partial  differentiation  with  respect  to  r  and  a 
dot  denotes  diffei entiation  with  respect  to  A^,  A^,  etc.  will  be  defined 
shortly. 

Let 


,  2 

W  *  )2n_i 

\y  / 

-"(.VO 

y 

4  1  -  ll)2" 

=  n(  1  -  J— ) 

\  2  / 

\  X/ 
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then 


°n  *  "0  - 

£)♦(*)*■-» 

(51-c) 

a  ..(4 

n  \y2 

(51-d) 

en  =[“2( 

y  x  /  \  n  xn  1  J 

(51-e) 

ii ;(®)/r  - 

“n  1.  I- n(l  -y2)  -  2y2  1  , 

r  \n  +  1 

n  'r*1'  ” 

4y2  I  **-(n+l)(l-y: 

2)j^ 

,x) 

+ 

- 

d  j~n(l-y2)  t  2y2~j  /  x\n-l 

1 

(52-a) 

«-► 

II 

Xe«  jap-twV* 

4y  f  nLin  +  l)(l-y2)r 

](!)' 

n  +  1 

- 

/  _r\n- 1 

V  x/ 

4 

-i 

(52-b) 

vn  - 1)(1  -y  )-* 

S^(8)(r,t)  - 

t  (*)■-* 

■„(£r2 

- 

nC„(7)n‘2tdn'n*2> 

©"! 

(52-c) 

S„9(8)(r.t) . 

■T  j-*„(n*2,(f)n+ 

nbn( 

.f)°  +  2 

4 

nCn(7)"'2-dn<”-2» 

(f)" 

(52-d) 

TjjS)(r,t) 

ne  (  ,  4  , 

a  (— )n+  b  (— ) 
2  |  n  \  x/  n  \  r/ 

n  +  2 

c„(7)n'2-d„(t)1 

(52-e) 
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The  functions  A^,  A°,  etc.  are  given  below  for  n  -  { 1 , 2 , 3,  .  .  .  ,N} . 


A  (co) 
n 


&JL.  .  J  4y>“  ),  2  4yV  Wl),  vvF(l)'«)  x  v 

2  2  4  2  A-0  2^/Fn  (^,X’>)Fn  V7’X,J, 

xv  x  y 

AV'  F^(u).x,y)F<1>(-a,x.y)J 

K  F'3l(M,x,y)F<1)(a,x,y)] 


n  n 
4  2  0 


2  2 


2  2 


.2  „  2 


“T"  -|V[F<1>,M,x,y,F<4\«,x.y).  F^W.ylF^.x.y^ 

X  \ 

A  ^r['f  )<^x-v»Fi3\f  >f  W,y>F<2)(f  .x.yj] 


X  y 
4  2r 


(53-a) 


A®(r ,<4  -  \  “Vyr  F^’tuj.x.r)  -  rjj 

X  J 

'-f  AV'  ^I^'toF.ylF^f.x.y)  Ff,u,,r,y)F«1\a  x,y)J 


~2  2y2  '  2  4y 2n2  ,  ^(1).  ,  „(3).  o» 


X 

2  4 


V 


^<lWr.y)F<4\f,x,>-).F<ni!»(«.r.y)lf\ax.,l] 

X  \ 


(53-b) 


A  •  (r  /J.V  2  ,22J?.r  iuj2  -  lF(l)ioe  2^F(3)(tt)  .1 

n'  ,UJ}  2  2  L  ^  2  /  n  \y,x,r/  ,2  n  \ ’ x ’  V J 

NX  X 

-A  Arv“2-  j2T-iFn\f'r'Vyfi,>«"’X.-', 

x  y 

-  -7-  -,r,v)I'"(“>)(uJ,x>v).  r(2)  ^,r,\  i  K<l)(u),x,  v)l 

•1  2  2  L  n  s>  ■  /  n  '  ’  ’•  n  \y  ■  i  n  '  J 


x  v 


TT 


X 

2.4 


2A^.  i^)^)(M  r,y)F(3Uxiy( 


MF^--^iVx,y,+  rf)(Mr,y)F(3)(WlX,y)J 


.  2  4  (  _  4  o  2 

-  17  y  )  4  8y  n  (n  - 1)  .2  2  2/  1  1  \1 

n  'W'  4  2  2  2  2  2  *  4y  n  u>  \-g  +—) 

wxy  (it  l  xy  V  /J 


(53-c) 


+  n 


r'4y4n2(n2- 1)  .2  2  2  1  1  \  4/x2  v2  \1 

~2*2  »"^+-2j+“(,— +  V)J' 

x  y  y  x  y  x  -1 


2r.  4  2,  2  „  .222 

x  I  4y  n  (n  -1)  4y  n  u) 

2  L  4  '2 

y  x  x 


n  ;(w»x»y)Fj1  \-,x,yj 

u ]F^)(W,x,y)Fji2)(-^,x,y) 
W4 ]FIi3)(W.x,y)Fji3)^,x,y) 
^  Fn4)( w>  x»  y) rjj4^* x ,  y) 

2/ 


„2  I-,  4  2,  2  „  .2  2  2 

JL-fe.  n  (n  - 1)  _  4*_  ry^  ^  ,  4]r(3) 

x2  L  y4  y2 

.4,2 


x  y 
2  2  2 


^[<1Wx,,)F;2)(a,x,y),lf)(Wlx,y)F|i1\^,x,y)J 

‘  i^K’Wx,y,..-<3,(M,iyj ,  F <3W, «<>(*„)] 

■  ^f[F®(l*x,,)F|14)(S|x,y)(  FW(to,x,y)Ff(a,x,yj] 

■  (u.x.ylF^a  x.,;+  r^Wx.yJF® ^,x,y)]j 


n <'•“>= -Vrr  fi^[F(„3Wx.r)  -  F,1»(«,x.r)J 

wx  y  (  *-  11  -1 

-  n2(u)2 x2  -  2y2n2)F<1\uJ,r,y)F<1)!^,xfy! 


(53-d) 


n  V7,A,-V 
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-  2y2n2F^({0,  r,y)F^^~,x,y) 

+  (A2-  2y2n2)F^(W,r,y)Fji2)fex,y) 


+  2y2F;i2>(a},r,y)FW(Mx,y)|  (53_e) 

B:<r,“)  '  '  iv  (^2-^2n>W(f,x,rj] 

+  2y2n2F^(Mr,y)F(l)(w,x>y) 

^yVrJ^rjjF^xj) 

-  2y2FJ12)(^r,y)Ff  (a),x,y) 

+  2y2Fji2)(-^,r,y)Fji4)(a;,x,y)|  (53-f) 

An/(r,w^  =T  r{^  'HT“[Fi2)(^x'r)-FJ,4)(^.x,r)] 


(53-f) 


.2  2  s  , 


(-^)^-^)Ff(.,r,y)F(1)(f,x,; 


.  2y2  (  2 

+  — hz —  ill 


^^.441  (3) 


r  V10  "2~  JF„  W r..v,F[|3)(^,x,y) 


(54-a) 

k",(r'“)  ^  ^[- 1“2-  5>v? 

2  2  2 

2y  n,  2  4y  n  \  J3),ti)  ,  m 

2  W  ~2  /  n  \  y  ’  r  ’  VJ  n  ( u)’ x  > v ) 
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4y 


2 

x  v 


r[F!i3)(y*’r*y)F!f)(w'x’y)4  Fn4)('y>riy)Fn1)(w’x’y)_ 


jy 

2 

x  y 


(54-b) 


f'[Fn3\rr'5')Fn4)<“’x'>'>+Fn4)if'r'  > J Fn3><M>x- 


B*>,u»  =  i2_n_[FW(oj,x,r)_FCT(ii),Xir)J 

iO  x  y  (  u  -1 

-  n2(a2x2-  2y2n2)Fjj3)( ou,  r, yJF^1  ^ x, y) 

-2y2n2Fji3)(W)r,y)F^3\f)x(: 
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4  2v2F^4V 


2y“F^^,r1y)Fji4^,x,y 


(54-c) 
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(5(i-l,) 


(56-d) 


2  y 

i'n'4>(‘°)'^  | F^2 ’ 4)( w>  +  F^3 •  4)( „)  +  n2 F^1  • 2 >( „)  +  n2 F^1  • 3 >( co) 

'  J¥_[7'Fn  ,(“'x'y)F!12)(rx'yj+ W*.y)J 
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-  [^IfWwIF'J^xj),  FW(«,x,y)F^(aiX,y^| 
Fn’3><“)  =  i  |I'i2'4'(“)  +  F<3’4»«o)H„2F<1’2),tt»,n2F<‘’3),u)  <Jl>‘e) 

‘  ^  [7^  »?W.y)F«(*,x.y)] 
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( 56— f) 
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if  4>«->  -  i  -  (*V  ±4-  -  »2)f^W  „2F<2'3><*, 


,  „2f^4>«4  -  ^i[A.F<3»(Wlx,y)Ff  (Mx.y 


F^V.x.y>f  (w.x,y)F'1,(i-®f,x,y' 

y~ 


4  Fji1)(co,x,y)Fj)4)(-^,x,y^ 

^4,4)(W)  -  ~  jnV2’4)(co)  ♦  n2Fji3’4)(to) 
2  2,-.  ... 


(56— i) 


"  J^[^2F!i4)(tt,,X’y)F!f\f,X,y)4  I<n)('y,x,y)Fn2)(a>,x,y)J 


(56— j) 


where 


f'J,(")-|[Fn)^'x-')Fn>(y’x'y) 

The  completes  the  analytical  solution  of  the  equations  of  the  three  dimensional 
theory  of  elasticity  for  the  stated  problem.  In  the  next  section  some  observations  of 
the  form  of  these  solutions  will  be  offered. 
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C.  Discussion  of  the  Elasticity  Solution 


The  most  conspicuous  attribute  of  the  solutions  given  by  (17),  (18)  and  (40) 
is  their  double  series  form.  The  sums  shown  in  (17)  and  (18)  are  over  all  angular 
harmonics  (n)  contained  in  the  load  representation.  Since  the  load,  by  earlier 
agreement,  contains  a  finite  number  of  harmonic  components  (N),  the  sums  (17) 
and  (18)  are  finite.  However,  from  (40)  it  is  observed  that  each  harmonic  com¬ 
ponent  or  mode  of  the  response  is  itself  an  infinite  series.  To  establish  the  con¬ 
vergence  of  these  series  a  brief  study  of  the  asymptotic  behavior  of  the  characteristic 

values  w, 5  and  the  characteristic  functions  W  V  .  etc.  for  large  values  of  the 
nj  nj  nj 

summation  index  j  will  nc.v  be  presented.  The  following  results  may  be  deduced 
from  the  asymptotic  behavior  of  the  cross  products  of  the  Bessel  functions  for 
large  values  of  &d  given  in  Appendix  II. 


D  (uj) - — —  sina/j 

°  JXY 

<57— a) 

Dj(cu)~  ~^y~  sin  3-u;  sin  -y 

(57— b) 

2 

D„(  od)  ^  ‘xy-  sin  ciu;  sin  -y  ;  n^  2 

(57-c) 

Since  the  characteristic  values  u;n-  are  the  roots  of  the  equation  Dn(od)  =  0  we 

conclude  from  (57)  that 

“oj~  "¥•  •  K*  •  • 

(58-a) 

Ktt  Ln 

anj~  \  ,  y  3  ;  (K,L)  --  1,2,3,.  .  . 

(58-b) 

For  any  given  value  of  n,  the  integers  K  and  L  are  linearly  related  to  j  for  all 

j  above  some  minimum  value  J(n).  The  formulae  (58)  will  be  very  useful  when 

calculating  the  characteristic  values  numerically.  Gazis  \  7l  also  obtained  the 

relations  (57)  in  his  stud\  of  hollow  cylinders.  Using  these  approximations  for 

:12 


k _ 


^nj  the  series  (40)  are  found  to  behave  asymptotically  for  large  j  as  follows. 
For  n  =  0,1,2, .  .  . 


W~  f  ^<r-Y,cos  f, 


~  i  Jf  ^Sin^(r-Y,oosi|l. 
SnJ<r)Qnj(,)  ~  O'2’'2)  S  Jf  ^stnf  ,r-Y)COSi|!t 


(59-a) 


(59-b) 


(59-c) 


For  n  =  1,2,3,.  .  . 


W[f  CM  ££<X-rH-l)K  i0os|I(r-Y)J  |  cos  fit 


V  ,(r)Q  /X  1  f x  Ltt 

nj  r3yY  r  l3  sinyLn  J^r  cos  a  (x-r) 


jjcosyL’Tcos 


T  (r)Q  ,(t)  ~  isnz 
nj'  '  nj'  '  ^2 


c°s-^(r-Y)J 

~K  <-l)K41  f  I 

r  K2 


cosyMt 


(59 -d) 


cos  ~^(r-Y) 


Kl  [  Y  8in  ^  (X_r)4  ^~x  sin  (r-Y)j|cos 


Tnj(r)Qnj(t)~  *2*2.  _1_  fl_  -|b  Ln 


(  x  cosyL-rsin 


*  tin  __J _ n  ,  Ln,v 

r  ^2  sin  y  L/t  I  Y  sin  a  (X-r) 

-tTsin (r-Y)J  cos  y  t 


(59 -e) 
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In  each  of  the  above  expressions  the  first  term  occurs  at  the  roots  —  and  the 
second  at  y  —■ . 

We  may  conclude  from  the  above  asymptotic  relations  that  all  the  series 

in  (40)  do  in  fact  converge  for  (Y  r  ^  X). 

In  addition  to  verifying  the  convergence  of  the  series  (40),  the  relations 

(59)  also  provide  a  useful  aid  to  understanding  the  nature  of  the  solutions.  For 

iJSll/r..Y-t  t) 

example,  it  is  noted  that  all  these  relations  are  either  of  the  form  e  a 
i  (r-Y  t-  yt) 

or  e  a  '  .  Both  represent  traveling  waves  having  in  the  first  case  a 

unit  phase  velocity  corresponding  to  the  dimensionless  dilatational  wave  3peed 

and  in  the  second  case  a  phase  velocity  y  corresponding  to  the  dimensionless 

shear  wave  speed.  Therefore,  we  may  interpret  (.>9-b,c)  as  yielding  discontinuous 

r  9 

dilatational  stress  waves  in  S^fr.t)  and  S  (r,t).  These' discontinuities  or  steps 
in  the  stress  are  the  result  of  the  dilatational  wave  produced  by  the  suddenly 
applied  load  being  reflected  between  the  boundaries  of  the  shell.  From  (59-e) 
we  see  that  the  modal  shear  stress  is  composed  of  both  dilatational  and  shear 
waves,  however  then  are  continuous  and  in  the  form  of  a  ramp  function  rather 
than  a  step  function.  This  is  to  be  expected  since  there  are  no  discontinuities 
in  the  shear  stress  introduced  at  the  boundaries  as  was  the  case  for  the  radial 
stress.  From  (59-a)  we  observe  that  the  radial  displacement  is  continuous  and 
dependent  primarily  on  the  dilatational  wave  while  the  circumferential  dis¬ 
placement  (59-d)  is  dependent  upon  both  the  dilatational  and  shear  waves. 

At  this  point  it  will  be  convenient  vo  consider  the  practical  computation  of 
the  series  (40)  in  a  specific  problem.  Since  these  series  converge  at  least  in 
the  manner  of  a  step  function,  they  may  be  terminated  after  summing  some 
finite  number  of  terms  to  obtain  an  approximation  of  the  desired  function. 

Just  as  in  the  load  representation,  the  Banczos  smoothing  process  may  be 
applied  to  these'  Unite  sums  to  increase  the  accuracy  of  the  approximation. 


This  process,  as  explained  in  Appendix  I,,  will  result  in  each  term  of  the  sum 

CO  u 

being  multiplied  hv  the  factor  — —  sin  — ' ■ —  where  j  K  is  the  last  term 

“K  ' 

retained  in  the  sum. 

It  is  interesting  to  note  that  the  form  of  the  solutions  (17),  (IB)  and  (40) 

is  exactly  the  same  as  would  be  obtained  bv  solving  the  problem  by  the  Williams 

or  mode  acceleration  technique  [  37l.  In  this  context  the  terms  etc. 

are  referred  to  as  the  "static"  modal  solutions,  the  u)  .  are  called  the  natural 

nj 

frequencies  of  the  system  and  W  . ,  V  etc.  the  eigenfunctions  of  the  system. 

M  ill  m 

In  fact,  for  the  suddenh  applied  load  the  "static"  solution  is  the  solution  of  the 
corresponding  static  problem  in  which  both  the  equations  and  boundary  conditions 
are  independent  of  time. 

This  concludes  the  anahtical  investigation  of  elasticity  theory  and  we 
will  now  proceed  to  examine  a  specific  example  for  the  comparison  of  the 
theories. 
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IV.  A  NUMERICAL  EXAMPLE 


In  order  to  obtain  a  better  insight  into  the  nature  of  the  elasticity  solution 
and  also  a  clear  comparison  with  the  shell  theories  a  specific  shell  geometry  and 
load  distribution  will  be  studied.  An  interesting  example  which  also  has  physical 
applications  is  the  suddenly  loaded  cylindrical  arch  shown  in  Figure  (1-B).  A 
radial  load  of  constant  intensity  PQ  is  suddenly  applied  to  the  outer  surface  at 
(-  8  <  0  <  0)  of  a  cylindrical  arch  whose  ends  at  0  =  +  77  are  free  to  move  hori¬ 
zontally  on  frictionless  rollers  but  restrained  from  moving  vertically.  The 
proper  boundary  conditions  at  the  ends  are  v  -  Tgr  0  at  0  =  ±  j-.  This  problem 
is  analogous  to  the  problem  of  a  complete  cylindrical  shell  of  identical  thick¬ 
ness  and  material  properties  subjected  to  the  same  loading  as  the  arch  plus 
the  symmetric  reflection  of  this  load  about  the  0  f  planes.  In  this  com- 

plete  shell,  because  of  the  symmetry  of  the  load,  v  t  q  0  at  0  -  ±-~ ,  The 

rw  l 

proper  load  distribution  function  for  this  example  is  therefore 

(1  ,  lei  ifi 

g(0)  {0  ,  3'  d  -  8 

(l  ,  '  -3  '  9  '  * 

The  Fourier  coefficients  of  this  function  are 


ao  ^ 

an  -  n  sin  n* 

an  0 
bn  -  0 


n  -  2 , 4 . 
n  1,3,5,  ..  . 
n  1,2,3,... 


The  thickness  ratio  used  in  this  example  is  t:  0.1.  The  only  material 

property  required  in  the  analysis  is  Poisson's  ratio  which  was  chosen  to  be 

1/  0.3  from  which  the  wave  speed  ratio y  -  0.5315  max  be  obtained. 
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The  solution  of  this  example  by  the  shell  theories  is  given  in  Appendix  III. 

Next  the  angular  extent  of  the  load  (£)  and  the  number  of  terms  (\)  us<  >1  n 
the  load  representation  must  be  chosen.  In  order  to  detect  all  the  peculiarities 
of  each  theory  and  to  emphasize  the  possible  differences  between  the  theories  it 
is  desirable  to  choose  $  as  small  as  possible.  However,  as  8  is  decreased  the 
number  of  terms  (N)  necessary  to  obtain  a  reasonable  approximation  of  the  load 
increases  very  rapidly.  Therefore  some  compromise  must  be  made  to  obtain  a 
sufficiently  concentrated  load  representable  by  a  reasonable  number  of  terms. 
After  studying  several  different  load  representations  it  was  decided  to  choose 
8  -  i  -  0. 1  with  N  100.  The  load  distribution  function  g^(0)  is  shown  in 
Figure  (!' ,  .ov  these  parameters. 

Having  specified  all  the  load  and  shell  parameters  we  may  now  proceed  to 
the  actual  computation  of  the  solution.  First  the  natural  frequencies  must  be 
found  for  each  theory.  For  the  shell  theories  this  is  a  straightforward  com¬ 
putation  which  involves  finding  the  roots  of  a  quadratic  polvnomial  for  th<  Mugge 
theory  and  a  cubic  polynomial  for  the  improved  theory ,  for  each  harmonic  n. 
However,  for  the  elasticity  theory  this  involves  finding  the  roots  of  the  trans¬ 
cendental  equation  r  0  for  each  harmonic  n.  If  approximate  values  lor 

these  roots  are  known  they  may  be  used  as  starting  values  for  a  first  order 
Newton-Raphson  iteration  scheme  to  obtain  the  natural  frequencies  to  any 
desired  accuracy.  The  asymptotic  relations  (5b)  may  be  used  in  this  pro¬ 
cedure  as  follows.  For  each  "alue  of  n,  Dn(u>)  is  plotted  for  o  until  it 
is  noted  that  the  roots  are  obeying  the  asymptotic  formalae  (5b).  The  starting 
values  for  the  lowest  roots  are  then  obtained  from  the  plot  and  the  starting 
values  for  the  remaining  roots  are  obtained  from  the  asymptotic  formulae. 

Table  (1)  lists  some  of  the  starting  values  for  the  first  few  harmonics  and 
indicates  when  the  asymptotic  relations  become  valid  by  using  the  appropriate 

asymptotic  formula  from  (58)  rather  than  its  numerical  \aluo. 
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STARTING  VALUES 


0 

1 

2 

3 

4 

5 

6 

0 

0.90 

TT 

2- 

3- 

4— 

&iT 

6— 

a 

a 

a 

a 

a 

a 

1 

1.27 

TT 

ya 

TT 

a 

CO 

p  Ip 

2~ 

“a 

.  n 
4ya 

2 

0.06 

2.01 

TT 

ya 

TT 

a 

o  n 

2ya 

3ya 

2- 

a 

3 

0.19 

2.84 

TT 

IT 

a 

2y— 

a 

2^ 

4 

0.37 

3.70 

TT 

v— 

•  a 

2 

CL 

2y— 

'  a 

_  TT 

3y- 
'  a 

2- 

Q. 

All  the  larger  values  of  n  follow  the  pattern  set  by  n  -  2, 3  and  4,  however  the 
value  of  j  at  which  the  pattern  becomes  evident  increases  as  n  increases.  This 
presents  no  difficulties  since  the  starting  values  for  n  6,8, 10,  ...  ,  100  may 
be  computed  from  a  difference  scheme  as  follows. 

wnj  “  “n-2,j  *  <“„-2,J  ' 

In  this  scheme  co  0  •  and  Hi  ,  .  are  the  correct  frequencies  obtained  from 
n-2,j  n-4,j 

previous  calculations  and  u)n-  is  the  starting  value  for  the  Newton-Raphson 
iteration.  The  difference  scheme  was  formulated  to  compute  only  the  even 
harmonics  since  the  odd  harmonics  are  not  used  in  the  present  example.  Using 
this  procedure  to  generate  the  starting  values,  the  natural  frequencies  were 
obtained  to  eight  significant  figures  with  usually  no  more  than  three  iterations. 

Note  that  the  j  -  0  root  for  n  =  1  seems  to  belong  with  the  j  -  1  group  of 
roots  for  :he  higher  harmonics  n  -  2,3,4.  This  occurs  because  the  lowest 
frequency  for  n  =  1  is  Hi  -  0  which  gives  rise  to  the  rigid  body  motions  already 
contained  in  and  equations  (48-a,b).  Thereioro  it  is  absent  from  this 
Table  and  it  is  not  a  root  of  D^(co)  -  0. 


39 


TABLE  II 

COMPARISON  OF  NATURAL  FREQUENCIES 


N 

j 

Elasticity 

Theory 

Improved 

Theory 

FIukro 

Theory 

0 

0 

.  90455 

. 90988 

. 90988 

0 

l 

5 1. 15019 

0 

•> 

G2. 85898 

0 

1 1 

94.25255 

0 

•i 

125.  GG727 

0 

:> 

1 .17.  08248 

0 

« 

1 88. 49792 

2 

0 

. 0G9G4 

.  00957 

. 00997 

2 

l 

2. 01489 

2.  01519 

2. 02001 

2 

2 

1G.  94201 

17.  9.9108 

2 

:s 

91.2G898 

o 

•i 

99. 8  moo 

2 

5 

GO. 98990 

2 

(» 

G2. 87G08 

•i 

0 

. 97019 

. 90970 

. 97942 

•i 

1 

9.  70998 

9.71459 

9.72048 

asp 

2 

17.28G1G 

17. 09410 

:t 

90.  90794 

♦t 

94 . 999 1 0 

1 

5 

90. 97718 

0 

G2.  98G77 

r, 

0 

.89090 

. 84940 

.90011 

(5 

1 

9.40990 

5.4794.9 

5.49799 

o 

o 

17.  8.9991 

18.  27920 

o 

*1 

.1 

90. 49G87 

o 

‘1 

.9.9.  01195 

6 

5 

90.  90979 

G 

G 

09. 10919 

10 

0 

2.  22094 

2.21974 

2.5G819 

10 

1 

8.97999 

9.  05090 

9.  08409 

to 

2 

19.41009 

19.95099 

10 

o 

«> 

29. 70484 

10 

■l 

90. 00987 

10 

5 

90.  945.92 

10 

G 

09.  79088 

20 

0 

6. 99599 

0. 89078 

10. 98795 

20 

1 

1 7 .  2 1 4 1 G 

18. 01759 

18. 10980 

20 

2 

24.97077 

20. 09098 

20 

:i 

29.  97554 

20 

4 

4 1 . 99574 

20 

r» 

50.  51445 

20 

G 

00, 19098 

no 

0 

22. 99459 

22. 08089 

45. 10705 

no 

1 

29. 02909 

44.  74  099 

05. 19109 

no 

2 

42.  05902 

49. 85975 

no 

•> 

•  J 

48. 12222 

no 

•1 

57.80819 

no 

5 

58. 92040 

no 

G 

79. 89090 

100 

0 

48. 01200 

48. 97410 

90.55017 

100 

1 

51.48581 

88.  01057 

200. 79590 

100 

2 

01.89290 

94.47102 

100 

;t 

79.  98704 

100 

4 

84. 79185 

100 

5 

99.98815 

100 

G 

102.  10479 
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A  partial  list  of  the  natural  frequencies  predicted  by  both  the  elasticity  theory 

and  thf  shell  theories  is  given  in  Table  (2).  Only  even  harmonics  are  shown  since 

the  odd  harmonics  are  not  used  in  the  example.  Note  that  the  lowest  frequency  is 

q  indicating  that  the  shell  offers  the  least  resistance  to  motion  in  the  n  --  2 

mode.  Since  the  n  =  2  mode  is  also  prominant  in  the  load  representation  we 

expect  this  to  be  the  dominant  mode  of  response  for  the  shell  in  this  example. 

The  period  of  the  response  in  this  mode  is  T  =  — -  =  90  dimensionless  units 

2 1  0 

of  time.  Next  observe  that  the  frequencies  predicted  by  the  Fliigge  theory  are 
consistently  greater  than  the  corresponding  frequencies  predicteci  bv  either  the 
improved  or  elasticity  theories  with  the  exception  of  n  -  0.  For  n  2  the  lowest 
frequency  predicted  by  the  Fliigge  theory  is  only  0.475%  greater  than  the  lowest 
frequency  predicted  by  the  elasticity  theory,  however,  this  difference  increases 
as  n  increases  so  that  at  n  -  100  it  is  87.59?.  Therefore,  the  Fliigge  theory  should 
satisfactorily  predict  the  response  caused  by  the  low  harmonic  components  of  the 
load.  This  characteristic  will  rapidly  deteriorate  as  n  increases.  On  the  other 
hand,  the  lowest  freouency  predicted  by  the  improved  theory  is  always  within  lr(' 
of  the  corresponding  frequency  predicted  by  the  elasticity  theon  for  the  range  of 


harmonics  covered  by  this  Table.  Therefore,  we  may  expect  the  improved  theory 
to  satisfactorily  predict  the  response  caused  In  all  the  harmonics  in  this  example. 
One  of  the  most  obvious  differences  between  the  theories  is  the  number  of 
frequencies  associated  with  each  theory  for  the  various  harmonics.  In  the 
elasticity  theory  a  complete  set  of  radial  eigenfunctions  and  associate  eigenvalues 


(frequencies)  is  required  to  represent  the  proocr  radial  variation  of  the  response. 
The  shell  theories  on  the  other  hand,  only  represent  the  gross  effects  of  these 
radial  variations  for  any  position  Q  on  the  shell.  Since  for  each  harmonic,  except 
n  -  0,  the  response  is  due  to  the  combined  effects  of  shear  and  dilatation,  the 
elasticity  theon  contains  two  sets  of  frequencies,  one  associated  with  the 
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dilatalional  effect  (oj  ,  —  )  and  one  associated  with  the  shear  effect 

I* -1 1  •  •  •  j  Cl  jii« 

Lj't  •• 

(w  .  y— —  ).  The  Fltigge  theory  approximates  only  the  lowest  frequenc\ 

nu  j  •  •  i j  J  j  i  •  i 

from  each  set  (G0nQ,  con^)  while  the  improved  theory  includes  also  the  second 
frequency  associated  with  the  shear  effect  (cdn0,  GOnj,  oon9). 

When  computing  the  response  predicted  by  the  elasticity  theory  for  each 
harmonic,  the  series  will  necessarily  be  truncated  after  summing  a  finite  numl>er 
of  terms.  From  (59)  we  see  that  at  the  boundaries  (r  -  X,  Y)  all  the  series  converge 
uniformly  for  all  t.  However,  for  Y<  r<  X  there  are  step  functions  or  dis¬ 
continuities  in  the  radial  and  hoop  stresses  which  periodicall..  recur  with  ponod 
t  =  2a  at  any  given  value  of  r.  To  accurately  represent  these  step  functions  the 
series  will  be  summed  with  the  Lanczos  smoothing  factor  up  to  and  including  the 
K  =  100  term.  This  corresponds  to  suming  Os  j  s  100  for  n  0  r\id  0  <  j  2hh 
for  u  5  2. 

One  further  interesting  feature  of  the  response  may  be  deduced  from  this 
Table  of  frequencies.  It  may  be  shown  F  34 j  that  for  improved  theory  for  large  n 
Wnl~  O.  -lOOn 

u>n2 °*  ^n 

go  „  ~  0. 929n 
nJ 

while  for  Fluggc  theory 

go  ,  -  0.020 In2 
nl 

go  0  ~  0.9035n 
nz 

If  we  simultaneously  examine  the  angular  and  time  dependence  of  the  response  for 

any  value  of  j  it  is  of  the  form  cos  go  t  cos  n  9  which  ma\  be  written  as 
go  .  •'  xv 

cos  n(9  f — ^-t).  In  the  improved  theory  lor  am  value  oi  i,  -  \'  ulvre 

n  ‘  n  i 

V.  is  a  constant  so  that  the  solution  is  a  sum  of  harmonic  components  the 

form  cos  n(9  r  Vt).  This  represents  a  wave  propagating  .around  tin  t  m  underence 

of  the  shell  with  phase  velocity  V.  Thus  it  is  possible  in  the  m.pr.rib  th«  .,r-  b>, 

observe  wave  phenomena  related  to  the  angular  coordinate  but  not  t  <  the  radial 
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coordinate.  In  the  Flixgge  theory  — - —  ~  .  0261n  so  that  the  harmonic  components 
of  the  response  are  now  in  the  form  cos  n(0  -t  .0261nt).  Since  the  phase  of  each 
component  varies  with  n,  the  high  frequency  components  will  be  out  of  phase  with 
each  other  and  a  traveling  wave  will  not  be  observed. 

The  elastiticy  theory  frequencies  listed  in  Table  (2)  were  checked,  when 
possible,  with  those  given  by  Armenakas,  Gazis  and  Herrmann  [19]  and  com¬ 
plete  agreement  was  found.  The  asymptotic  relations  [58]  were  derived  earlier 
by  Gazis  [7]  and  [ll]  and  further  discussion  of  free  vibration  characteristics 
may  be  found  in  references  [6]  through  [19]. 

We  now  have  all  the  information  necessary  to  quantitatively  compute  the 
response  of  the  shell  predicted  by  each  of  the  three  theories.  The  results  of 
this  computation  are  presented  in  Figures  (3)  through  (15).  The  comparison  of 
the  theories  is  shown  for  the  radial  and  circumferential  displacements  at  the 
median  surface  of  the  shell  (r  =  1)  and  for  the  hoop  stress  at  the  inner  and  outer 
surfaces  (r  =  0.95,  1.05).  The  radial  and  shear  stresses  predicted  by  elasticity 
theory  at  tne  median  surface  (r  =  1)  have  also  been  computed. 

A  comparison  of  the  static  solutions  is  shown  in  Figures  (3),  (4),  (5)  and 
(6).  Note  the  excellent  agreement  between  the  improved  and  elasticity  theories. 

At  their  maximum  values  the  improved  and  c’  ist  city  theories  differ  by  less 
than  half  a  percent  in  these  results,  howev-i  ?  ':s  Flugge  and  elasticity  theories 
differ  by  five  to  six  percent.  The  maximum  stat  radial  displacement  occurs  at 
0-0.  The  maximum  static  tangential  displace:  4  occurs  at  0  -  0.70775  and  the 
maximum  static  hoop  stress  occurs  at  0  0,  r  0.95.  At  any  given  value  of  r 

and  0,  the  maximum  value  of  the  dynamic  response  is  expected  to  be  twice  the 
static  value  (see  [4],  pp.  181-182).  Therefore,  the  displacements  and  the  hoop 
stress  predicted  by  the  three  theories  in  the  dynamic  case  will  be  compared 
the  values  of  r  and  0  given  above. 
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The  static  shear  and  radial  stresses  predicted  In  the  elasticity  theory  at 
the  median  surface  (  r  1.00)  are  shown  in  Figures  (7)  and  (s).  The  maximum 
stress  occurs  in  both  cases  near  the  edge  of  the  load  at  0  0. 1 1715  as  expected. 

The  maximum  magnitude  of  the  hoop  stress  is  seen  to  be  twenty  times  greater 
than  the  maximum  shear  stress  and  forty  times  greater  than  the  maximum  radial 
stress.  This  provides  justification  for  the  usual  assumption  made  is  shell  theory: 
the  radial  stress  may  be  neglected  compared  to  the  hoop  stress.  We  will  now 
proceed  to  examine  the  dynamic  response  of  the  shell. 

Figures  (9)  through  (12)  depict  the  initial  response  of  the  shell.  In  this 
very  early  stage  of  the  response  the  various  stress  waves  may  be  observed  as 
they  propagate  through  the  shell  and  are  reflected  from  its  boundaries.  The  time 
history  of  the  radial  stress  at  the  median  surface  directly  beneath  the  load  is 
shown  in  Figure  (9).  The  events  depicted  on  this  graph  may  be  explained  as 
follows.  At  t  -  0  the  radial  stress  at  the  outer  surface  is  discontinuous!  y 
changed  from  zero  to  minus  one.  This  discontinuity  in  the  radial  stress  pro¬ 
pagates  as  a  compression  wave  into  the  shell  with  unit  dimensionless  velocity 
(dilatational  wave  speed).  At  t  =  .05  we  observe  this  compression  wave  as 
it  passes  the  median  surface.  At  t  -  t  it  encounters  the  inner  surface  of  the 
shell  and  since  this  surface  is  stress  free  it  is  reflected  as  a  tensile  wave  which 
is  observed  as  it  passes  the  median  surface  at  t  -  0. 15.  At  t  2 a  this 

tensile  wave  encounters  the  outer  boundary  from  which  it  reflects  as  a  com¬ 
pression  wave.  This  phenomenon  is  repeated  periodical!!  with  period  T  2x. 

The  time  historv  of  the  stress  at  the  median  surface  becomes  more  complicated 
as  each  wave  passes  and  adds  its  effect  to  those  of  the  previous  waves.  This  is 
the  reason  tor  the  changing  form  of  ti.e  response  curve  in  Figure  (9). 

Figure  (10)  shows  the  shear  stress  predicted  by  elasticity  theorv  at  the 
median  surface  of  the  shell  and  at  the  edge  of  the  applied  load.  As  predicted 
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earlier,  there  are  no  discontinuities  in  the  shear  stress,  however,  the  various 
waves  may  still  he  detected  since  they  cause  discontinuities  in  the  slope  of  the 
curve.  At  t  =  0.  05  we  observe  the  dilatational  wave  and  then  at  t  ~~  = 

0.935  we  observe  the  shear  wave  as  both  pass  the  median  surface.  Then  again 
at  t  =  0. 15  and  t  *0. 187  the  dilatational  and  shear  waves  reflected  from  the 
inner  surtace  are  observed. 

Figure  (11)  depicts  the  time  history  of  the  hoop  stress  at  the  median 
surface  directly  beneath  the  load  at  A  0  A  comparison  of  the  three  theories 
is  presented  on  this  plot.  In  the  elasticity  theory  the  dilatational  wave  is 
observed  as  it  is  reflected  between  the  shell  boundaries.  As  noted  earlier 
this  wave  cannot  be  predicted  by  the  shell  theories;  however,  they  do  accurately 
characterize  the  average  value  of  the  stress. 

Figure  (12)  shows  the  hoop  stress  predicted  by  the  three  theories  at  the 
inner  surface  of  the  shell  at  0  -  .  Because  of  the  wave  character  of  the 

elasticity  theory,  no  response  is  observed  until  the  dilatational  wave  originating 
at  the  edge  of  the  load  reaches  this  location  on  the  shell.  This  time  interval  is 
approximately  1-^-0  1.47.  Similarly  the  improved  theory  contains  waves 

traveling  around  the  shell  with  phase  velocities  V  O.  llMi,  V9  0.878  and 
r  0.929.  Therefore,  there  is  no  response  observed  until  the  fastest  wave 

^-8 

V  passes  this  point  on  the  shell.  This  occurs  at  t  —  =  1.58.  As  indicated 

.  3 

previously  the  Fliigge  theory  response  is  not  entirely  composed  of  traveling  waves 

and  therefore  it  predicts  an  immediate  response  at  every  point  on  the  shell. 

Figures  (13)  through  (15) show  one  full  period  of  the  response  predicted 
by  the  three  theories.  The  radial  and  circumferential  displacements  of  the 
median  surface  of  the  shell  are  shown  in  Figures  (13)  and  (11).  Note  the 
excellent  agreement  between  the  improved  and  elasticilx  theories.  In  both  eases 
they  differ  by  less  than  half  a  percent  in  the  vieinih  o!  the  maximum,  The 


Flugge  theory  however,  differs  from  the  elasticity  theory  by  approximately  seven 
percent  in  this  vicinity.  The  static  solution  predicted  by  the  elasticity  theory  is 
aiso  shown  on  these  graphs  and,  as  expected,  it  is  approximately  one  half  of  the 
maximum  value  of  the  dynamic  response.  Also,  the  period  of  the  motion  is  seen 
to  be  T  =  90  as  predicted. 

Figure  (15)  shows  the  hoop  stress  at  the  inner  surface  directly  beneath  the 
load  predicted  by  the  three  theories.  Again  the  agreement  between  the  elasticity 
and  improved  theories  is  excellent  both  in  magnitude  and  form.  Since  the  response 
predicted  by  the  Flugge  theory  is  slightly  out  of  phase  with  the  other  two  theories 
(the  Flugge  theory  predicts  a  slightly  faster  response)  there  are  large  differences 
between  the  theories  at  any  given  instant  of  time,  however,  the  maximum  value 
predicted  by  the  Flugge  theory  differs  by  only  nine  percent  from  the  maximum 
predicted  by  the  elasticity  theory. 

The  radial  and  shear  stresses  were  computed  for  one  complete  period  of 
the  response  and  their  maximum  magnitude  was  found  to  be,  as  in  the  static 
case,  less  than  four  percent  of  the  maximum  hoop  stress.  Thus  the  assumption 
made  in  shell  theory;  i.e. ,  the  radial  stress  is  negligible  compared  to  the  hoop 
stress,  is  also  valid  in  the  dynamic  case. 

This  concludes  the  study  of  the  specific  example.  The  principal  findings  of 
this  investigation  will  now  be  summarized  in  the  conclusion. 
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FIGURE  3.  STATIC  RADIAL  DISPLACEMENT  VS.  POLAR  ANGLE 


120 


HGURE  6.  STATIC  HOOP  STRESS  \/S.  POLAR  ANGLE 
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FIGURE 
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FIGURE  11.  INITIAL  RESPONSE:  HOOP  STRESS  VS.  TIME 


FIGURE  12  INITIAL  RESPONSE  =  HOOP  STRESS  VS.  TIME 


FIGURE  15.  RADIAL  DISPLACEMENT  VS.  TIME 
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FIGURE  15.  HOOP  STRESS  VS.  TIME 


V.  CONCLUSIONS 


i 

i 


» 


The  principal  findings  of  this  investigation  will  now  be  summarized.  There 
are  three  convenient  stages  of  the  response  for  which  the  shell  and  elasticity 
theories  have  been  compared.  First,  there  is  an  initial  response  in  which  the 
dilatational  and  shear  waves  transmit  the  effects  of  tl  e  loading  to  the  various 
points  in  the  shell  and  the  radial  variation  of  the  response  begins  to  develop. 
Neither  of  the  shell  theories  can  accurately  describe  the  details  of  the  response 
during  this  early  period  and  the  elasticity  theory  must  be  employed  to  observe 
this  phenomenon.  The  period  of  the  wave  phenomenon  occurring  during  this 
initial  response  is  T  2  i  which  is  the  time  required  lor  a  dilatational  wave  to 


travel  from  the  outer  surface  to  the  inner  surface  and  then  back  again  to  the 


The  second  stage  of  the  response  consists  of  the  effects  of  the  load  being 

transmitted  around  the  shell  in  the  circumlerential  direction.  The  characteristic 

periods  of  the  response  produced  bv  the  dilatational  and  shear  waves  in  this 

‘> 

stage  are  nppioxi match  T  2  and  T  ~.  As  shown,  the  improved  theory 
does  provide  a  good  approximation  of  this  .aspect  of  the  response  with  the' 
periods  being  slightly  larger  than  those  predicted  In  elastieitv  theory.  However, 
the  Flugge  theory  cannot  characterize  this  aspect  ol  the  response  because  of 
the  nature  ol  its  frequency  spectrum. 

The  third  stage  of  the  response  contains  the  maximum  displacements  and 


stresses  occurring  in  the  shell  tor  all  time.  One  complete  period  ol  the  response 

♦) 

occurred  in  the  interval  T  — - — -  whereof  was  the  lowest  frequence  predicted 

ui,  2,o  i.i 

2,o 

by  the  various  throw  s.  As  explained  earlier  tin-  n  2  mode  is  the  principal 


harmonic  component  of  the  response.  Since  the  Flugge  theory  very  accurately 
characterizes  this  mode  of  response,  the  results  predicted  by  the  Flugge  theory 
are  surprisingly  good.  The  inaccuracy  in  the  Flugge  theory  characterization  of 
the  higher  modes  contributed  to  a  total  error  of  less  than  10^  in  the  maximum 
hoop  stress  and  even  less  in  the  maximum  displacements.  However,  the  improved 
theory  is  still  far  superior  to  the  Flugge  theory  because  of  its  accurate  characterization 
of  the  higher  modes.  Thus  the  error  in  the  improved  theory  was  less  than  1/2$  in 
both  the  maximum  stress  and  displacement.  Also  it  was  seen  that  the  I  liigge  theory 
always  underestimated  the  maximum  stress  whereas  the  improved  theory  provided 
a  slightly  conservative  estimate. 

If  the  applied  load  is  modified  in  such  a  way  that  the  relative  importance  of 
the  higher  modes  is  decreased  the  Flugge  theory  may  be  expected  to  yield  an  even 
better  approximation  of  the  overall  response.  This  would  occur  if  the  load  were 
applied  at  a  finite  rate,  if  the  angular  extent  of  the  load  were  increased  or  if  the 
load  were  distributed  continuously  over  the  shell  surface.  On  the  other  hand  if  the 
load  becomes  more  concentrated  the  error  in  the  Flugge  theory  response  will 
increase.  However,  since  the  relative  difficulty  of  both  shell  theories  is  about 
equal  compared  to  elasticity  theory  and  since  the  improved  theorv  yields  a  closer 
approximation  to  the  actual  response,  the  improved  theory  must  be  judged,  io  be 
superior  to  Flugge*  shell  theorv. 

One  turther  lactor  to  be  considered  when  comparing  the  shell  theories  to  the 
elastici.,  theorv  is  the  computation  time  necessan  to  obtain  the  solution.;  All 
computations  lor  the  example*  presented  were  programmed  in  Fortran  IV  on  a 
CDC  computer.  The  computing  time  required  lor  each  shell  theorv  was 
approximated  ten  minutes,  however  the  computing  time  required  for  the 
elastmitv  theorv  was  approximated  lour  hours.  Klasticity  theory  required  a 
greater  .amount  <•!  time  because  ot  the  number  of  terms  summed  in  the  series  lor 


01 


each  harmonic  to  obtain  the  initial  response.  If  the  initial  response  is  not  desired, 
considerably  fewer  terms  may  be  used  in  these  series.  However,  the  computing 
time  will  still  greatly  exceed  that  required  by  the  shell  theories  because  the  Bessel 
functions  involved  in  the  solution  require  more  computing  time  than  the  algebraic 
functions  involved  in  the  shell  theory  solution.  Thus  if  the  initial  response  is  not 
desired,  this  cons.deration  represents  an  additional  advantage  of  improved  shell 
theory. 

In  conclusion,  it  appears  that  improved  shell  theory  is  the  superior  choice  of 
the  three  theories  in  terms  of  the  trade-off  between  accuracy  and  analytical  com¬ 
plexity  if  the  details  of  the  initial  response  are  not  required.  If  the  initial  response 
is  desired  elasticity  theory  must  be  used.  We  also  note  that  there  are  several 
shell  theories  which  are  similar  to  Flu’gge's  theory.  Thus  it  may  be  expected  that 
the  present  findings  with  respect  to  Fliigge's  theory  are  also  applicable  to  the  shell 
theories  of  Love,  Donnell,  Vlasov,  Sanders,  etc. 
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APPENDIX  I:;  LANCZOS'  SMOOTHING  TECHNIQUE 


The  accuracy  of  the  approximation  of  a  function  over  some  interval  by  a 

truncated  Fourier  series  may  usually  be  increased  by  applying  the  operation  of 

local  smoothing  to  the  function  being  represented.  This  procedure  along  with  some 

illustrative  examples  are  presented  by  Lanczos  in  his  book  on  Fourier  series 

([38],  pp.  61-75).  The  crux  of  the  method  is  the  replacement  of  the  original 

function  by  01  which  is  locally  smooth.  For  example,  suppose  we  wish  to 

represent  the  function  f(X)  defined  over  the  interval  -L<  X<  L  by  a  Fourier 

th 

series  which  has  been  truncated  after  the  K  term  in  the  series.  Then 

%  r 

f  (X)  =  ~  +  )  (a.  cos  ul.X  -  b.  sinoo.X) 

is.  &  /-,  )  J  J  J 

j;  1 

where  aQ,  a^  and  b.  are  the  Fourier  coefficients  of  f(X)  and  ^  Consider  a 

new  function  f(X),  derived  from  f(X)  as  follows. 

jt_ 

f“K 

f(X)31^  J  ^Y)dY 

n 

“k 

The  following  conclusions  may  be  obtained  through  the  application  of  the  mean  value 

theorem  to  the  above  integral.  If  f  is  continuous  at  X  ,  then  lim  [  lim  f(X)l  ~  f(X  ). 

0  X-*X  K-»  »  0 

o 

If  f  is  discontinuous  at  X  ,  then  lim  [  lim  7(X)]  — — where 

0  X-»  X  K-»  "  2 

o 

f '  lim  f(X(  ■  “  )  and  f“  lim  f(X  -  o).  For  large  but  finite  values  of  K,  this 
-t  -*  o  '  ->  o 

smoothing  procedure  has  very  little  effect  on  f(XQ)  if  f  is  continuous  over  the 

interval  (X  -  —  c  X  X  •  — — ).  However,  if  f  has  a  discontinuit\  at  X  , 

°  0  0 
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the  discontinuity  will  be  smoothed  into  a  rapid  but  continuous  change.  The  truncated 
Fourier  series  approximation  of  this  smoothed  function  may  be  obtained  by  sub¬ 
stituting  into  the  integral  definition  to  obtain. 


a0  ~  ooK  u;.rr 

fK(X)  "  T~  T  L  7J7T'  sin  ~£~  (ai  C0S  WiX  ‘  bi  sin  WiX) 
j  -  1  J  K  J  J  J 


The  accuracy  of  th.s  approximation  of  f(X)  is  better  than  the  truncated  Fourier 
series  at  all  points  except  those  in  the  immediate  neighborhood  of  a  discontinuity. 
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APPENDIX  n:  CROSS  PRODUCTS  OF  BESSEL  FUNCTIONS 

The  properties  of  the  cross  products  of  the  Bessel  functions  defined  in 
equations  (28)  and  (42)  are  listed  below. 

For  0<X<  ®,0<  Y<  K  =  (1,2, 3,4)  and  n  =  0,1,2,  .  .  . 

(1)  ^(P.X.Y)  are  entire  functions  of  P 

(2)  FjiK)(-P,X,Y)  =  FjiK)(P,X,Y) 

(3)  Fj)K)(ioo,X,Y)  =  --|FjiK)(a),X,Y)  ;  (K, n)  4  (4,0) 
F<04)(iw,X,Y)=jFj>4)(w,X,Y) 


where  i  and  a)  is  real. 


lim 
P-*  c 

F<1) 
>  0 

<P,X,Y)  = 

1 

2n 

lim 
P-*  c 

F<2) 
»  n 

(P,X,Y)  - 

1  1 

"  2  1 

(i)"] 

lim 
P-»  c 

Ft'3) 
,  n 

(P.X,  Y)  = 

1 

2 

s¥r<m 

lim 
P-*  c 

F<4> 

,  n 

(P*X,  Y)  = 

-f[ 

i  _  i 

►<!>< 

3 

1 

t  3t 

lim 

P  -»  r 

p(4) 

,  0 

(P.X,  Y)  ^ 

■H 

T - 1 

(5)  As  P-»  ■»  the  functions  behave  asymptotically  as  follows. 

P;1)(P,X,Y)-  ^  sinhP(X-Y) 

v/XY 

F^  (P,X,  Y,~  -  yicosh  P(X-Y) 
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F^P.X.Y)-/!  cosh  P(X-Y) 
^(P.X^-VXY  P  sinh  P(X-Y)  ;  0 

^(P|X|Y)~  — - ^sinhP(X-Y) 

TxF  1 


(6) 

(a)  For  n  -  0 


Ti  Fo1)«""X*Y>  “  i  [Ff  Wx,Y)  *  F£V,X,Y)] 


F<f>(U,X,Y)  =  u)  Y [xF^4' (10  ,X, 7'  -  YF<,V,X,Y>] 

JZ  F™ (0),X,Y)  =uiX [yfW (uj , X, Y)  -  XF^ (u ,X, Y)] 

&  fW(1o,X,Y)  =  f  [fF^.X.Y)  +  |  F<3>(a,,X,Y)  ♦  2fW(«,X,Y)] 


(b)  For  n  =  1,2,3,  .  .  . 


k  Fa)<“,'X-Y»  '5  [f!,2)<“.x.y)  ♦  ifV.X.Y)] 

fJ,2)(".X,Y)  =I[fW(W,X,Y)  4  (n2-1U2Y2)F^>(a),X,Y)J 
-F®(u,X,Y).  ^[f!,4)<U.X,Y)h  (»2-»¥)fW(ui,X,Y)] 

fW(U,X,y)  ,  ^[laV  -  „2)  F|f>(U,X,Y)  <  (wV-nVJV.X.Y)] 


d 

doc 


J 

aoo 


a  co 


(?) 

(a)  For  n  -  0 


^jF<1\u,X,Y>-  jF^do.X.Y) 
Tx  F^loi.X.Y)  w2Y  fW(UiX,Y) 
"7Y  Fq2 ’ (c*J ,X, Y)  -  ^rf'lw.X.Y) 
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I 


gpyF<3)(W,X,Y)=w2XFjl)(<ofX,Y) 

(b)  For  n  -  1,2,3,  .  .  . 

Ac  ^(W.X.Y)*  1  F^co.X.Y) 
^F<k2>(w,X,Y)=  ^4>(W,X,Y) 

Ac  Fn3>("-X*' Y>  =  "  4-(w2X2  -  n2)  F(1) 


Ac  Fi4)<w*X* Y)  =  -  £  (a2*2  -  n2)F2(w 

^yF^W.X.Y)*  y  F^2>(u),x, Y) 

Fy  Fn3)(w,X,Y)  -  y  F^(jo,X, Y) 


L 


(W.X.Y) 

,X,Y) 


67 


APPENDIX  III:  SHELL  THEORY  SOLUTIONS 


The  response  of  a  cylindrical  shell  in  plane  strain  produced  by  an  arbitrary, 
radially  directed  load  P(0,t)  is  given  in  [39!,  p.  22  and  [ 34J,  p.  297,  lor  both  the 
Fliigge  and  improved  shell  theories.  In  terms  of  the  nondimensionalization  used  in 
this  study  these  solutions  may  be  rewritten  in  the  following  form.  Eor  the  Flu'gge 
Theory:' 


w(0 ,  t) 

v(0,t) 


N  2 

~~  Hr  E  (t)  t  >  '  L 

a-f  L2  o 

n  1  i  1 


A.  F.  (0,t)  1 
n  in  in1  ' J 


N  2 

1 

arr  L  <_ 
n  -  1  i  -  1 


L  B.  G.  (0,t) 
n  in  in'  ' 


V'9’1' 

(i-p; 

N  2 

.  "  V 


N  L  (A.  <  nil.  >'  *  n“(r-l>A.  )F,  (0,t)  j 

_  iv  in  m  '  ;ir  in' 

n  1  i  1 


For  the  Improved  Theory: 
w(0,t)  -  -  [i  Eo(t) 


N  3 
V 

nil  1 


L  A.  E.  10, t) 
n  in  in'  ' 


N 


v(0,t) 


.\ 


,  l,  B. 


n  1  i  I 

V -Oh 


B.  G.  (0,t)  I 
n  in  in'  J 


.  JL_  E  (t) 

r  2  9  1  o'1' 

(l-l >)  L 


L  (A.  nB  n(r-i)C.  )E.  <9,t)  ! 
_  n'  ip  in  '  in  m  '  j 


•  1  i  I 


lih 


* 


where 


F  (t)  -  f  f  P(  ■  ,T  )sin  w0(t-t  )d ;  dT 

°  ^0  0J  -TT 

Fin(9,t)  =  f 

G.  (0,t)  =  f  f  P(-  ,T)sin  n(8  -  ■lOsin'jdj  (t- T)d>  dr 


in  o' 


T  N  .  nrr 

L  -  - Sin  -rr 

n  nT  N 


•J 1  -  2 1/  ^  _  "J  o 

<0,_  -  ■  /T“.~  d.  i  -  /-,  ..v  il 


in  (1-y)  in  ’  "o  (1-1/ )  *  o 
For  the  example  discussed  in  Section  IV  P(9,t)  -  g(0)H(t) 


where 


(1  ,  O<0<  0 

g(-0)  ;  g(0)  =  <  0  *  fi<  8<  h-0 

(  1  ,  ’T-0<  0  <’T  , 


Therefore,  for  this  example 

r  1  -  COS  CO  t  -j 

roW'<— S - J 


ou„ 


F.fil0 ,  t)  -  ^  sin  n0  cos  n 0 


1  -  cos  UJ.  t 

m  o  i  /• 

J  ;  n  -  2,4,(>,  .  .  . 


uO- 

in 


4  f1  "  cos  ^  1  1 

Gin(0,t)  -  -  sin  n0  sin  n0  j - - - J  ;•  n  -  3,4,6,  .  . 


CO. 

in 


F.  G.  -  0  ;•  n  -  1,3,5. 
;n  m 


The  muial  coefficients  A.  ,  B.  ,  C.  and  the  natural  frequencies  Cl. 

in  in’  in  in 


soeicded  with  each  theory  are  given  belc’v.  For  the  Flegge  theory. 

p-  p 

0.  /t  .  a  ■  ,  a>„  ./ -r-d-xj 


(it) 


as- 


!  /3  '  X  2  < 

A.  =  JL  ( _ n  _  n  \  .  -1 

ln  X  \  2  T/’  A»-  = 

n  n 


2n  X 


n 


1-X  2 

_ n  n 

2  \ 

n 


Bln  \  X  *  B2n  = "  Bin 
n  n  in 


3X2 

n 


V  1  +  ^  +  ^(n2-!)2  ,  X„  =  J  1  - 

For  the  improved  theory: 

■ 

A‘" ' 

Bin  d~  La5n  ^in  “  Wln^ 

Cl  '  <i~  ^a6nglc_  a7n*fn] 

i2  i  r  -2 


in 


in 

o 


din  *  C^i,n'^in]  C*i+2,n*4]  M0D  <3> 


a2  4i  r  , _ 

^11  '  9  £  1  “  •/l-4g  J 

’  M  01 


2  =9,  _ 

'  ~2~  1  -t  '1-4^^  J 

^31  '  0 


.2  1  r  *2(i-l)-r 

£.  —  <>■  i  _  9 v  n  '  \ 

yin  9  hi)  ^  —A  COS - - - — —  i 

in  J  -ii  i_  n  \  3  J 


,  n  *•-  2 


xn  *I-3s.-,„ 


arc  cos  i-j  !  -  -  <>  ~i„  if,  „  ,)r.  2 

"  I  -  sHn  2  -nrjL1 


'i.'l 


c  -  r.r  f 


2  2 
a  n 


V 


g2n  =  %“<2+K  ><1+-l2->  +  K‘  (1+  -%-) 


2  4 
a  n 


,2W,  ^3a 


2  2 


\4a' 


g3n  “  nT  (1  +2K  X1  f  12  ^  +  K  n  (1  +  -  V) 


+  K2  a+ii-)2(1+4-) 


a2K2 


o-  -  -  ••  /i  ,  a  .  2.  2  ,,2 
g4n  ”  12  (1  +  12  *  n  <h  “ 


’5n 


glng3n 


’2n 

2„2 
a  K 


glng4n 


36n 


g, 


12 


2  +  K 


L=4  +  V  + 

Z  4 


2n 


2  i  _  JL. 

144 


K 


5  +  5 _ 

x- m  r2  (:  zj  ) 


'1  +  2K 


K4  = 


1  -f  K 
K2 


2  ,2 

a_ 

12 


In 


z  >  ( 

n 


,,-2  2 


OQT  2  /  ‘  \1 

i?  K  v1  12  yj 

2  . 


a 


4)  *«2(’  4)( 


i  + 


■2- ) 
4  J 


a  =.  — —  (  i  _  — _ 

on  12  12 


2  2 


I-  4  6  Z 

a4n‘n|?iir  ,K  ll4Tr)J 

a5„ '  4r  [>  *  2«2  1 4  <>  •  x2ij 


‘.in  K  ^  'f2- >('“-!) 


l7„  ”[K\ 


2  2  2  2 

"■h-;  '  TTV2  ‘  TT/1] 
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K  =• 


A  -  V  2 

J  2  K  where  K  is  the  Mindlin  shear  coefficient.  K 


example  of  Section  IV. 


=  0. 86  in  the 
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